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ABSTRACT 

Aims. To show the importance of high-spatial resolution observations of H n regions when compared with observations obtained with 
larger apertures such as ISO, we present mid-infrared spectra of two Magellanic Cloud H n regions, N88 A and N160 A. 
Methods. We obtained mid-infrared (8 - 13 yum), long-slit spectra with TIMMI2 on the ESO 3.6 m telescope. These are combined 
with archival spectra obtained with the Infrared Spectrograph (IRS) onboard the Spitzer Space Telescope, and are compared with 
the low-spatial resolution ISO-SWS data. An inventory of the spectra in terms of atomic fine-structure lines and molecular bands is 
presented. 

Results. Concerning N88 A, an isolated Hn region with no adjacent infrared sources, the observations indicate that the line fluxes 
observed by ISO-SWS and Spitzer-IRS come exclusively from the compact H n region of about 3" in diameter. This is not the case for 
N160 A, which has a more complex morphology. We have spectroscopically isolated for the first time the individual contributions of 
the three components of N 160 A, two high-excitation blobs, Al and A2, and the young stellar object N160 A-IR. In addition, extended 
[Siv] emission is observed with TIMMI2 and is most likely associated with the central star cluster located between Al and A2. We 
show the value of these high-spatial resolution data in determining source characteristics, such as the degree of ionization of each 
high-excitation blob or the bolometric luminosity of the YSO. This luminosity (2 x 10 5 L e ) is about one order of magnitude lower 
than previously estimated. For each high-excitation blob, we also determine the electron density and the elemental abundances of Ne, 
S, and Ar. 

Key words, infrared: ISM - ISM: lines and bands - ISM: dust, extinction - ISM: H n region - ISM: individual objects: N88A - ISM: 
individual objects: N160A 



1. Introduction 

High-Excitation Blobs (HEBs) are compact H n regions that con- 
stitut e a rare class of ionized n ebulae in the Magellanic Clouds 
(e.g. Heydari-Mala yeri et all I20071 and references therein). 
These objects are tightly linked to the early stages of mas- 
sive star formation, when the stars disrupt their natal molecu- 
lar clouds. They are considered to be the first step in the evo- 
lution from an ultracompact H n region towards more extended 
structures. Their main characteristics are: high excitation, small 
size (~ 2 pc), high density, and higher extinction than the typical 
Hn regions in the Magellanic Clouds. So far, only 6 members 
have been discovered in the Large Magellanic Cloud (LMC) and 
3 in the Small Magellanic Cloud (SM C). Most of these HEBs 
have a counterpart observed i n radio (Indebetouw et "afl l2004t 
iMartm-Hernandez et aHl2005l) . 

Infrared spectroscopy is one of the best ways to deter- 
mine the characteristics of H n reg ions. The overall mid -infrared 
(MIR) spectra of H n regions (e.g iPeeters et al.L l2002b) are dom- 
inated by atomic fine-structure lines that originate in the ion- 
ized gas, a strong rising continuum due to thermal emission of 
dust and strong, broad emission features at 3.3, 6.2, 7.7, 8.6, 
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1 1.2 and 12.7 /vm, attributed to emis sion by Polycyclic A romatic 
Hydrocarbon (PAH) molecules (e.g Allamandol a et allll989h . 

The comparison of ISO-SWS spectra of H n regions in the 
Magellanic Clouds (which also include a few HEBs) and in our 
Galaxy shows a clear segregation in PAH ratios between the 
sources of different type of envi ronment, i.e. from the Milky 
Way to the LMC and to the SMC (I Vermeil et aUl2002bl) . This is 
most certainly due to the different environment encountered in 
the Magellanic Clouds, in particular, its low metal content. This 
study, nevertheless, is based upon observations made with the 
large ISO aperture (between 14" and 33"), where other contri- 
butions besides that of the H n region itself are likely to exist. 

MIR spectra of Hn regions give also access to various el- 
ements, e.g. Ar, S and Ne, in two different ionization stages. 
Ratios of fine-structure lines such as [ Ar ra]/[ Ar n] 9.0/7.0, 
[Srv]/[Sm] 10.5/18.7 and [Ne ra]/[Ne n] 15.5/12.8 fxm probe the 
ionizing stellar spectrum between 27.6 and 41 eV and are widely 
used to constrain the properties of the ionizing stars in H n re- 
gions. Again, the use of large apertures is limited since the con- 
tributions from neighboring sources cannot be disentangled. 

High-spatial resolution observations are thus essential to 
provide accurate measurements of the line (and PAH) fluxes 
emitted by the object itself. In this sense, the study of Hn re- 
gions in the Magellanic Clouds, and specially of HEBs due to 
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Table 1. Journal of the TIMMI2 observations. 



Source 


RA (* m s ) Dec (° ' ") 


Date 




Standard 


Flux (mjy) 




(J2000.0) 






star 


Nl N10.4 


N88A 


01 24 07.9 -73 09 04 


2003-11- 


14 


HD4815 


100 ± 40 (a) 


N160A1 


05 39 43.3 -69 38 54 


2003-11- 


17 


HD39523 


1330 ± 270 (b) 


N160A2 


05 39 45.9 -69 38 39 


2003-11- 


■17 


HD39523 


1630 + 330 (b) 


N160A-IR 


05 39 43.9 -69 38 34 


2003-11- 


14 


HD39523 


< 40 (c) 



(a) The image was calibrated using the spectrum since we did not image the standard star in this filter. Slit losses of up to 50% are considered in 
the uncertainty. <b) The uncertainty (~ 20%) accounts for variations in the transparency of the atmosphere and the absolute flux calibration error of 
the standard star. <c) Upper limits are defined as the flux of a point source with a peak flux three times the image noise. 



Table 2. Journal of the Spitzer observations. The coordinates of 
the positions observed with Spitzer-IRS are given. 



Table 3. Journal of the ISO observations. The coordinates of the 
positions observed with ISO-SWS are given. 



Source 



RA (* "' ! ) Dec (° ' ") 
(J2000.0) 



AORKEY <a 



IRAC 






N88 A 


4383488 




N160A 


12546816 




IRS 


SH/LH 


SL/LL 


N88A 1 24 08.10 -73 09 05.00 


4384000 




N160A2 5 39 46.10 -69 38 36.00 


12606208 


8242944 




12902556 


8568832 




13048576 


8579328 




13162752 


11812352 




13358336 


12003584 




13529088 


12130048 




15672064 


15672064 




16295936 


16295936 




16702976 


16702976 


MIPS 






N160A 


14381312 





1 AORKEY which uniquely identifies each Spitzer observation. 



their compactness, might serve as a stepping stone to understand 
H ii regions in other galaxies, where their emission is confined 
into one beam. H n regions in the Magellanic Clouds are close 
enough that we can still resolve their individual components. 
Hence, we can disentangle their contribution to the overall spec- 
trum and assess how much the properties of the H n region de- 
pend on aperture size. This is actually difficult to do for Hn re- 
gions in our Galaxy where the surrounding ionized gas, because 
of its low density, is usually too diffuse and extended to be well 
sampled. 

In this paper, we study in detail the spectral morphology of 
two well distinct regions, N88A and N160A, at high-spatial 
resolution using TIMMI2 and the Infrared Spectrograph (IRS) 
onboard the Spitzer Space Telescope. These spectra are subse- 
quently compared with findings (and conclusions) based upon 
large aperture data obtained with ISO-SWS. N88A, located 
in the SMC, is a compact HEB whic h is among the brightest 
objec ts in the Magellanic Clouds (e.g. iHevdari-Malaveri et all 
1999). N160A, located in the LMC, is a particularly inter- 
esting region since it harbors tw o HEBs named Al and A2 
dHevdari-M alaveri & Testor, 1986) and a massive Youn g Stellar 
Object (YSO; Epch tein et alTll984HHenning et al.lll998l) . called 



Source 


RA (* '" s ) 


Dec (° ' ") 


Date 


TDT (a) AOT 




(J2000.0) 






N88A 


01 24 08.15 


-73 09 03.2 


1997-09-08 


66300110 01 (b) 


N160A2 


05 39 46.12 


-69 38 36.6 


1997-08-05 62804104 01 (b) 


N160A-IR 05 39 43.74 


-69 38 30.4 


1997-08-01 


62401303 01 (c) 


N160A1 


05 39 43.27 


-69 38 51.4 


1996-07-10 


23600925 02 



<a) Target Dedicated Time (TDT) which uniquely identifies each ISO 
observation. (b) Speed 2. (c) Speed 3. 



N160 A-IR. For N160 A, we have spectroscopically isolated for 
the first time the individual contributions of the two HEBs 
and the YSO . The images in the near- and mid-infrared of 
iHenning et alj (1 19981) clearly separate this YSO from the neigh- 
boring sources, but the spectrum published by these authors, ob- 
tained with ISO, is dominated by emission from the surrounding 
ionized gas. 

The core of this paper is structured as follows. Section [2] 
presents the observations and data reduction. Section[3]describes 
the general morphology of N88 A and N160 A. Results are pre- 
sented in Sect. [4] which includes a description of the MIR spec- 
tra, a comparison between observations obtained with different 
apertures, an estimate of the physical conditions of the ionized 
gas (in terms of electron density and elemental abundances) and, 
finally, the modeling of the Spectral Energy Distribution (SED) 
of the YSO in N160 A. The main conclusions are summarized in 
SectE 



2. Observations and data reduction 

2.1. TIMMI2 

The A^-band imaging and spectroscopic observations were made 
with the Thermal Infrared MultiMode Instrument (TIMMI2) at 
the ESO 3.6 m telescope. A journal of these observations is pre- 
sented in Table Q] The camera uses a 320 x 240 pixel Raytheon 
Si:As array. In imaging, we selected the ff.'2 pixel scale provid- 
ing a total field of view of 64" x48". For long-slit spectroscopy, 
we used the 10 pm low-resolution grism which ranges from 7.5 
to 13.9 pm and has a spectral resolving power A/AA ~ 160. 
The slit used was l'.'2x70", with a pixel scale of 0V45. It was 
oriented along P.A.= 0° covering the brightest nebular spots of 
N88 A (cf. Fig.CQ and N160A (cf. Fig.©. The three slits cov- 
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ering N160 A are denominated, from west to east, si, s2 and s3. 
They were selected to cover, respectively, N160 Al, N160 A-IR 
and N160 A2. The infrared seeing during our observations var- 
ied between l'/6 and 1'.'9, measured from the full-width-at-half- 
maximum (FWHM) of the standard star spectra. 

In order to correct for background emission from the sky, the 
observations were performed using a standard chopping/nodding 
technique. For N88 A, we used chop and nod amplitudes of 15" 
for imaging, and of 10" for spectroscopy. To image N160 A, we 
used an amplitude of 30"; for spectroscopy, we used 20" in the 
case of slit s3, and 30" in the case of slits si and s2. 

The data processing included the removal of bad frames and 
the co-addition of all chopping and nodding pairs. This proce- 
dure is applied to both the target and the calibration standard 
star. We propagated the uncertainty for each pixel, dominated by 
variations of the sky transparency, along each step of the pro- 
cessing. 

The A^-band imaging of the standard stars was performed in 
the same filter as the target and used for photometric flux con- 
version from photon count rates (ADU/s) into astronomical units 
(Jy). The MIR flux densities of the targets in the resulting images 
are given in Table[TJ N88 A was observed in the Nl filter (red and 
blue cuts are 7.79 and 9.46 //m, respectively), and N160 A in the 
N10.4 filter (9.46 and 11.21yum, respectively ). Astrometry was 
performed by comparing with radio maps (cf. Indebe touw et all 
2004; iMartm-H ernande z et all 120051) and by matching the posi- 
tions of stars in common in our images and in the 2MASS survey 
catalogue. The final (contour) images are shown in the left pan- 
els of Figs.[TJand|2] 

The spectra of the standard stars were used for atmospheric 
corrections and for flux calibration. To minimize residuals of the 
sky line correction, we observed the standard star and the tar- 
get at similar airmasses. In addition, we applied a flux correction 
provided by the TIMMI2 instrument team that takes into account 
the differential extinction between the target and the standard 
star. The synthetic calibrated spectra for these standard stars are 
given by Coh en etaO ([199 9). We note that, close to the blue and 
red cut-off frequencies of the grating, the terrestrial background 
varies strongly with ambient conditions and atmospheric correc- 
tions are more difficult. The same holds for wavelength regions 
of strong terrestrial lines near 9.58 pm, 1 1.73 /mi and 12.55 pm. 
Finally, wavelength calibration is straight forward since a table 
with the pixel-to-wavelength correspondence is provided by the 
TIMMI2 team. 

Spectra were extracted using apertures of 1.5 times 
the FWHM measured along the spatial direction, 
where FWHM(N88A)= 3'.'0, FWHM(N160 Al) = 3'.'1, 
FWHM(N160A2)= 5'.'3 and FWHM(N 160 A-IR) = 273. The 
resulting spectra are shown in Figs. [3] Hand |5] 

2.2. Spitzer 

N88 A and N160A were also observed with the I RAC instru- 
ment dFazio et all 120041) . the M IPS instrument dRieke et all 



120041) . and the IRS instrume nt (iHouck et all 120041) onboard 
the Spitzer Space Telescope dWerner et all 120041) . From the 



Spitzer archive, we obtained the IRAC pbcd data associated 
with pipeline version S14, the MIPS pbcd data associated with 
pipeline version S16, and the IRS Short-High (SH), Long-High 
(LH), Short-Low (SL) and Long-Low (LL) bed data associated 
with pipeline version S15 (see Table [2j. Some observations, al- 
though available, were not considered here. The LL spectrum of 
N160 A2, for instance, is not used since the LL wavelength re- 
gion is covered already by the SH/LH modes, which also have 



a higher spatial resolution. In addition, the LL (1 st order) spec- 
trum of N 160 A-IR and the 24 pm MIPS map of N 160 A were 
not considered due to saturation problems. 

IRAC images of N88 A and N160 A are shown in the right 
panels of Figs.[TJand[2] The IRS SH and LH apertures centered 
on N88 A (only nod 1) are indicated in the right panel of Fig.Q] 
Here, only the spectrum of nod 1 is used because of the better 
overlap between the different orders. N160 A2 was used as a cal- 
ibrator for IRS and therefore, it was observed multiple times us- 
ing the same pointing (see Table |2). The orientation of each slit 
depends on the observing date/time. The right panel of Fig. [2] 
for instance, shows the slit of one specific SH and LH observa- 
tion centered on N160 A2, i.e. AORKEY = 15672064. The other 
slits are rotated with respect to this one. IRS SL and LL spec- 
tra of N160 A2 were also obtained through multiple slits. This 
and the large slit size allowed us to create spectral maps using 
CUBISM -see e.g. Fig. [7] where a [S iv] 10.5 pm line flux map 
has been created using this software- and, consequently, to ex- 
tract a spectrum for each individual component. These extraction 
apertures are shown in the left panel of Fig. [2] Full-slit aperture 
extraction was performed for the SH and LH observations. The 
aperture sizes are listed in Table |4] 

Additional processing was applied to the IRS bed data. 
Rogue pixels and bad pixels identified in the bmask were re- 
placed with IRSCLEAN. This program obtains the new flux val- 
ues by deriving profiles for a group of rows above and below 
the bad pixel (i.e. neighboring wavelengths), predicting the pro- 
file for the bad pixel's row, and scaling that profile by the good 
pixels in the bad pixel's row. Residual bad data were manually 
removed from the obtained spectra. Note that background sub- 
traction has not been applied. The resulting SH and LH spectra 
have a resolving power of 600. The SL and LL spectra, on the 
other hand, have a resolving power ranging from 60 to 128. They 
are shown in Figs. [3] [4] and [5] 

Given that N160A2 is extended, different slit orientations 
result in flux level variations up to 20% in the resulting individ- 
ual SH and LH spectra (because of the full-slit aperture extrac- 
tion). For the final line fluxes reported in Table [4] we consider 
the mean of the fluxes in all positions. Uncertainties take into 
account these variations from slit to slit. 

In addition, in all SH 2-dimensional (2d) bed images of 
N160 A2, there are between 1 and 4 pixels in the [Ne m] 15.5 pm 
fine-structure line that only have 1 useable (i.e. non-saturated) 
plane in the raw data making it impossible to obtain a reliable 
flux for the considered pixel. These (saturated) pixels are cor- 
rected for by IRSCLEAN, which operates on the 2d bed images. 
However, the accuracy of the corrected flux decreases with the 
increasing number of bad pixels. Nevertheless, when comparing 
the [Nem] 15.5//m line flux of N160A2 for all SH observations, 
no correlation is found with the amount of pixels that have 1 sin- 
gle useable plane, and all line fluxes fall within +10cr of the 
mean value, similar to the range in fluxes found for the other 
fine-structure lines. This suggests that the variations we see in 
the [Ne m] line flux are most likely due to the different slit orien- 
tations (as mentioned previously). The [Nem] 15.5/vmline flux 
should anyway be considered with caution. 

2.3. ISO-SWS 

N88A and N160A we re observed with t he Sh ort Wavelength 
Spectrometer (SWS, Ide Graauw et all 1 19961) onboard the 
Infrared Space Observatory (ISO, iKessler et all 119961) in the 
AOT 01 and AOT 02 modes with a resolving power (A/AA) rang- 
ing from 350 to 1500 (see Table [3]). The AOT 01 mode covers 
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Fig. 1. (Left ) Contours of the Nl (A c — 8.70 fim) image of N88 A. Contour levels are 2.5, 4 and 6 times the image rms (cr = 0.05 
mJy/pixel). The ISO-SWS (band 2C, between 7.0/mi and 12.1 pun) aperture (dotted lines) and the TIMMI2 slit (dashed lines) are 
indicated. (Right) Spitzer-IRAC image of N88 A at 8.0 pm. The Spitzer-IRS SH and LH apertures (only nod 1) are plotted by solid 
lines. The dotted box outlines the field-of-view covered by the left panel. North is up and East is to the left. The (0", 0") position 
in both panels corresponds to a(J2000)=l h 24 m 7 s .920 and £(J2000)=-73°9'4 S .07 [see electronic edition of the Journal for a color 
version of this figure]. 



the full SWS wavelength range from 2.3 to 45/im, while the 
AOT 02 mode consists of separate line scans. The AOT 01 data 
were processed with IA 3 , the SWS Interactive Analysis package 
using calibration files and procedures equivalent with pipeline 
version 10.1 . A detai led account of the reduction can be found in 
Peeters et alj d2002bl) . For the AOT 02 observations , we re fer to 
Vermeii et al.l d200'2al) and lVermeii & van der Hulstl (120021) . The 
ISO spectra of N160A2 and the YSO N160 A-IR are plotted in 
Fig- El Note that the ISO-SWS observation centered on Al was 
made in AOT 02 mode and consists only of separate line scans; 
hence, this spectrum is not shown. 



3. Source characteristics 

3.1. N88A 

N88 A is among the brightest objects in th e Magellanic C louds 
and is part of the large Hn complex N88 (lHenizdll956t other 
designation is DEM 161. lDavies et al.lfl976h . located east of the 
main body of the SMC. It is a comp act (about 3" in diame- 
ter) a nd high excitation H n region (cf. iHevdari- Malaveri et all 
I1999I) . unusual among SMC Hn regions because it contains 
much dust fe.g. lTestor & PakullL[l985tlKurt et allll999h . The in- 
terstellar reddening towards N88 A, as derived from the Balmer 
decrement in the visib l e, is on average Ay ~ 1 .5 mag (cf. 
IHevdari-Malaveri et all Il999t iTestor et all 120051) . unexpect- 
edly high for an Hn region in the metal poor SMC. Hence, 
IHevdari-Malaveri et all d 19991) describe N88 A as a HEB, a very 
young Hn region which is hatching from its natal molecular 
cloud and is heavily affected by absorbing dust associated with 
the cloud. 



N88 A has a complex and inhomogeneous morphology (cf. 
IHevd ari-Mala veri et all 1 1999b . An absorption lane crosses the 
nebula from north to south and here the extinction rises to more 
than 3.5 mag. West of this dusty structure lies the brightest part 
of N88 A, a small core of 073 in diameter. A sharp edge is vis- 
ible on the north-west and several lower excitation arc-shaped 
features and filaments seem to emerge from the main body of 
the so urce. Its radio continuum emission (cf . Ilndebetouw et all 
2004) agrees fairly well with the global optical morphology, al- 
though it appears to be slightly offset to the west of the Ha emis- 
sion. These authors note, however, that this offset could be due 
to the ~ 1 " pointing uncerta inty of HST 

IHevdari-Malaveri et al.l (l999) also show that the bright Ha 
core of N88 A contains at least two stars (#1 and #2), and a third 
fainter one is located to the east, marginally detected and thus 
more embedded. T hese authors suggest that this faint star might 
be the exciting star. Tes ter et all ((2003) derived a spectral type of 
06-08 V for the ionizing star. This spectral type classification is 
strengthened by th e nebular near-infrared spectrum obtained by 
ITestor et all (|2005), where the ratio of the 2.11 /mi He i to Bry 
is consistent with an ionizing source of spectral type earlier than 
07 V (see table 3 of lHanson etail 120021) . 

Figure Q] shows the TIMMI2 (Nl filter) and Spitzer-IRAC 
(8 /an) images of N88 A. On top, we show the TIMMI2 slit, and 
the Spitzer-IRS and ISO-SWS larger apertures. 

3.2. N160A 

N160 dHenizel Il956t other d esignations are DEM 2 84, 
iDavies et alj 1 19761 and MC76, iMcGee & Miltonl 1 19661) is 
a giant Hn complex lying about 30' south of 30Doradus. 
The complex is formed by two main regions ~3'.5 apart 
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Fig. 2. (Left) Contours of the TIMMI2 N10.4 (A c = 10.38 ^m) image of N160 A. Contour levels are 2.5, 4, 7 and 10 times the 
image rms (cr = 0.21 mJy/pixel). The red stars (cf. Heydari-Malaye riet alll2002h show the positions of the brightest stars ionizing 
N160 Al (star #22) and N160 A2 (stars #100 and #105), and of the brightest stars grouped in the central cluster (stars #24, #35 and 
#58). The ISO-SWS (band 2C, between 7.0 pm and 12. 1 pm) apertures (dotted lines), the three TIMMI2 slits (dashed lines), and the 
Spitzer-IRS SL (blue solid lines) and LL (orange dashed lines) apertures centered on each component are indicated. The TIMMI2 
zones defined in Fig.|6](l, 2 and 3) are also indicated (black solid lines). (Right) Spitzer-IRAC image of N160 A at 8.0 pm. One of 
the Spitzer-IRS SH and LH apertures (solid lines, which correspond to AORKEY= 15672064) is plotted as a reference. The dotted 
box outlines the field-of-view covered by the left panel. North is up and East is to the left. The (0", 0") position corresponds to 
a-(J2000)=5 h 39 m 43 s .291 and £(J2000)=-69 o 38'54:28 [see electronic edition of the Journal for a color version of this figure]. 



dHevdari-M alaveri & Testorl H986). The brighter western region 
comprises components A (NGC 2080) and D (NGC 2077), while 
the eastern region is formed by components B (NGC 2085) 
and C ( IC2145). N160 is linke d with the OB star association 
LH 103 dLucke&Hodgd 119701) . 

N160A is the brightest component of N160. It is a partic- 
ularly interesting region sinc e it harbors two HEBs (cal l ed Al 
and A2, and discovered by iHevdari-Malaveri & Testorl [1986) 
and a young stellar object (named N160A-IR). This YSO 
was discovered using near- and mid-infrared photometry by 
Epchtein et al. ( 198 fj[) and has been further observed in the in- 
frared bv lHenning et al.l d!998l) . A water maser has been found 
at the interface between the molecular cloud and the H n region 
providing strong evidence that feedb ack from massive st ars is 
triggering subsequent star formation dOliveira et all 120061) . 

HST observations of N160 A dHevdari-Malaveri et al.Ll2002h 
reveal a bright nebula of about 35" by 25". The compact Hn 
region Al, which is the brightest part of N160 A, shows a tiny 
cavity of some 2'.'3 across which might be the result of the strong 
wind of a bright star found towards the center of the blob. This 
star must be a massive O type since it has carved the cavity 
and produces the highest [Om]/Hj6 ratio in the whole region. 
The other HEB, A2, is 3" in diameter, has a patchy appearance 
marked by the presence of several thin absorption lanes and con- 
tains about dozen rather faint stars. Three bright cluster stars are 
also identified within the nebula along with the two blobs Al 
and A2. These three stars seem to contribute to the ionization of 
N160A. The visual magnitude of the brightest member of this 
central cluster agrees with an 06.5V star. 



The Al and A2 blobs lie in very dust rich areas of the re- 
gion. In the case of Al, dust is mainly located behind the south- 
ern border of the cavity where the mean value of the HafH/3 
Balmer decrement gives an extinction Ay = 1.3 mag. The north- 
ern compact H n region A2 is more affected by dust, with a mean 
extinction Ay = 1.5 mag (cf. lHevdari-Malaveri et al.L 12002). 

The radio continuum emission observed towards N160 A re- 
veals an extended structure coincident with the optical nebula 
where the tw o HEBs Al and A2 stand out as two strong ra- 
dio p eaks (cf. Ilndebetouw et"aQ 12004; Mar tin-Hernandez et all 
2005). As expected from its youth, the YSO does not have asso- 
ciated radio continuum emission. 

Figure|2]shows the TIMMI2 (N10.4 filter) and Spitzer-IRAC 
(8^m) images of N160A. The YSO is not registered by the 
TIMMI2 image but it appears as a very bright source in the IRAC 
one. The larger field of view of the Spitzer-IRAC image also 
registers the western N160D region. The TIMMI2 slits, and the 
Spitzer-IRS and ISO-SWS apertures are indicated. 



4. Results 

4.1. The MIR spectrum 

We present here the high-spatial resolution spectroscopy data 
obtained with TIMMI2 and Spitzer, which will be compared 
with ISO observations. 

A variety of fine-structure lines and broad-band dust fea- 
tures fall within the ground-based A^-band spectroscopic range 
(8-13 pm) of TIMMI2. The most relevant lines are [Arm] at 
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Fig. 3. Spteer-IRS SH spectrum of N88A. The TIMMI2 spec- 
trum is plotted in grey. The most important atomic lines and PAH 
bands are labeled. Note that this Spitzer spectrum is not scaled 
down. 



9.0pm, [Siv] at 10.5pm and [Nen] at 12.8pm. These lines re- 
quire hard radiation with energies between -21 and 35 eV, and 
the most likely explanation for their excitation mechanism is 
photoionization by hot stars. Spitzer-IRS and ISO-SWS cover a 
wider wavelength range, including lines such as [Ar n] at 7 pm, 
[Ne in] at 15.5 and 36 pm, and [S m] at 18.7 and 33.5 pm. 




12 14 16 

Wavelength (/inn) 

Fig. 4. Spitzer-IRS SH spectrum of N160A2. The TIMMI2 
spectrum is plotted in grey. Note that the Spitzer-IRS spectrum 
is not scaled down. 



In terms of molecular features, the MIR spectra of many H n 
regions are dominated by the well known PAH emission features 
at 3.3, 6.0, 6.2, 7.7, 8.6, 11.0, 11.2 and 12.7pm. Two silicate 
bands centered around 9.7 pm and 18 pm can also be present. 

The TIMMI2, 5p;'fzer-IRS and ISO-SWS line fluxes, and the 
upper limits of the non-detected lines, are listed in Table [4] The 
PAH band fluxes are reported in Table [5] 
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Spitzer-IRS SL: N160A-IR 




ISO-SWS x 0.1: N160A-IR 

TIMMI2 i 




(e) 



6 8 10 12 14 

Wavelength (fim) 

Fig. 5. MIR spectra of N160 Al, A2 and the YSO N160 A-IR. 
Panels a, b and d show the Spitzer-IRS SL spectra in black and 
the TIMMI2 spectra in grey. These Spitzer-IRS spectra are not 
scaled down. Panels c and e compare the TIMMI2 spectra of 
N160A2 and N160A-IR (in grey) with the (scaled down) ISO- 
SWS AOT 01 spectra centered at these positions (in black). 



4.1.1. N88A 

Figure [3] shows the TIMMI2 and Spitzer-IRS SH spectra of 
N88 A. The spectral appearance of the ISO-SWS spectrum 
matches that of the Spitzer-IRS spectrum and therefore, we do 
not show it. 

N88 A is resolved by TIMMI2 with a FWHM along the spa- 
tial direction of 3'.'0 (see Sect.|2~TT), higher than the TIMMI2 slit 
width of 1'.'2. Hence, slit losses are expected. Assuming that the 
source is Gaussian and that the slit is perfectly centered on the 
source peak, the slit might be registering about 40% of the to- 
tal emission of N88 A. Assuming this slit loss, the total [S iv] 
line flux associated with N88 A would be ~ 85 x 10~ 20 W cm -2 , 
in good agreement with the fluxes measured by Spitzer-IRS and 



ISO-SWS. The same holds for the dust continuum flux. This in- 
dicates that, in spite of its larger aperture, the fluxes observed 
by Spitzer-IRS, and even ISO-SWS, come exclusively from the 
compact source N88 A. This is expected since N88 A is an iso- 
lated object with no known, adjacent infrared sources. Thus, it 
can be said that Spitzer observations do not spatially resolve 
N88A. 

It is interesting to note the lack of PAH emission features in 
the TIMMI2 spectrum. The Spitzer-IRS spectrum, on the con- 
trary, clearly exhibits a 1 1 .2 pm PAH band. Other PAH bands at 
e.g. 11.0 and 12.8 pm, and the PAH emission between 15-20 pm, 



N. L. Martin-Hernandez et al.: Mid-infrared study of the H n regions N88 A and N160 A 



7 



are however not present in this Spitzer spectrum. In addition, hy- 
drogen recombination lines at 11.3 and 12.4 /mi are detected. 

4.1.2. N160A 

The HEBs A1 and A2 As it can be seen in the top three panels 
of Fig.[5](and in Fig.|4|, the spectra of Al and A2 are the typical 
spectra found in H n regions, characterized, as discussed above, 
by emission from fine-structure lines and PAHs on top of a dust 
continuum. 

Both sources are resolved by TIMMI2, with a FWHM along 
the spatial direction (see Sect. |2.1t higher than the TIMMI2 slit 
width. Again, slit losses are expected, in particular for A2 with a 
FWHM of 5'.'3. In principle, one could expect that the TIMMI2 
slit is registering about 40% of the total emission in the case of 
Al, and only 25% in the case of A2, assuming that the sources 
are Gaussian. However, the TIMMI2 and the Spitzer-IRS spec- 
tra of these two sources match very well (see panels a and b 
of Fig. [5]) suggesting that slit losses are not very significant, at 
least not for the dust continuum emission. A comparison with 
the scaled down ISO-SWS spectrum of N160 A2 is shown in the 
middel panel of Fig. [5] The large difference in the flux level of 
the continuum between the TIMMI2 and the ISO-SWS spectra 
reflects the effect of the much larger ISO-SWS aperture, and the 
extensiveness and complexity of the source and its surroundings. 

Aperture differences are also evident in the line fluxes mea- 
sured by each instrument (see Table @). In the case of Al, for 
instance, the [S iv] line fluxes measured by TIMMI2 and Spitzer 
agree with each other, but they only account for 20% of the [S iv] 
line flux measured by ISO. For A2, the [S iv] line flux measured 
by Spitzer is in between those obtained with TIMMI2 and ISO, 
which clearly reflects the effect of the extended emission asso- 
ciated to this whole area. In particular, the [S iv] line flux mea- 
sured by TIMMI2 accounts for approximately 40% of the [S iv] 
line flux measured by Spitzer. This percentage is higher than the 
expected 25% (see above) and hence, we are confident that the 
smaller aperture of Spitzer-IRS (i.e. 4'.'7xll") is indeed tracing 
most of the ionized gas associated with A2. 

When looking at the [Ne m]/[Ne n] 15.2/12.8 pm line ratios 
measured for N160A (see again Table 0), aperture effects are 
also evident. Note, for instance, how this ratio (and also the 
ground-based [Siv]/[Nen] 10.5/12.8/vm ratio) changes depend- 
ing on the position and/or aperture. Basically, these ratios tend 
to decrease with aperture size. This has a clear impact if one 
wants to use this ratio as a diagnostic for the effective tempera- 
ture of the central ionizing star or stars. The [Ne m]/[Ne n] line 
ratio is highly sens itive to the eff e ctive t emperature of the ion - 
izing star (see e.g. iMokiem etall l2004t iMorisset etail [2004), 
and an incorrect line ratio could lead to a wrong interpretation 
in terms of the stellar content. When sampling a single H n re- 
gion, the use of pencil beam spectra might give a -sometimes- 
wrong impression of the ionization characteristics of the central 
star. Essentially, the low ionization stage of an ion is undersam- 
pled for lines of sight towards the center of the region, and the 
other way around for lines of sight away from the center. In this 
sense, it is much better to use observations that encompass the 
whole Hn region. However, N160A is not a single Hn region. 
It harbors two HEBs (each containing at least one massive O 
star) and the main ionizing stars (from which the hottest could 
be an 06.5 V) appear to be part of a central cluster which ex- 
tends over several arcseconds. This cluster is located between 
Al and A2 (see the left panel of Fig.|2|. The ionization structure 
of N 160 A is, therefore, highly complex. The variations we see 
in the [Ne m]/[Ne n] (and [S iv]/[Nen]) line ratio, with the high- 



est value found close to the HEBs, suggest that the two HEBs 
are ionized by very hot stars, while the extended emission is 
ionized by cooler stars that are likely located in the central clus- 
ter. These two HEBs could be ionized, in fact, by hot 04-05 V 
stars, as estimated from their 4. 8 GHz radio continuum emission 
dMartfn-Hernandez et aill2005l) . 

Finally, N160A1 exhibits PAH emission bands at 6.0, 6.2, 
6.8, 7.7, 8.6, 1 1.2 and 12.7 pm (see panel a of Fig . [5) . It has a re- 
markably strong 6.0 pm band which starts already near 5.83 pm. 
The 6.2/11.2 and 7.7/11.2 band ratios (see Table |5]l are slightly 
lower than those of N160 A2, but similar to those found towards 
LMC Hn regions ( I Vermeil et all l2002bl) . In addition, the band 
profiles of the main PAH bands belong to class A and are, there- 
fore, con sistent with the typica l PAH characteristics found in H n 
regions dPeeters et all |2002a). The top of the 7.7 pm complex, 
however, is fairly sharp and emission seems to be lacking from 
the 7.8 pm component. This could possibly be due to the less 
reliable fluxes at the order edges and to mismatches between the 
different orders of Spitzer-IRS SL (one order overlap falls ex- 
actly at the position of the 7.7 pm complex). 

N160A2, on the other hand, exhibits PAH emission bands 
at 5.75, 6.0, 6.2, 7.7, 8.6, 11.0, 11.2 and 12.7 /mi (see panel b 
of Fig. |5]). The profiles of the main PAH bands (at 6.2, 7.7, 8.6 
and 11.2 pm) belong to class A. As for the fine-structure lines, 
aperture effects are apparent. The PAH band fluxes measured by 
Spitzer are considerably smaller than the ISO ones. In particular, 
they only account for about 24% and 15% of the, respectively, 
7.7 and 1 1.2 pm PAH fluxes measured by ISO. This percentage 
is roughly similar to the corresponding aperture ratio, suggesting 
that the PAH emission is very homogeneous. In addition, there 
exists a significant difference in the 7.7/(11.0+11.2) band ratio 
between the Spitzer and the ISO data. Given that the 7.7 and 1 1 .2 
pm emission is dominated by PAH cations and neutrals respec- 
tively, this suggests a higher neutral to cation ratio in the larger 
ISO-SWS aperture. In general, the 7.7/(11.0+11 2) band ratio 
is sim ilar to those found in LMC H n regions by Verme iTet alj 
(2002bJ), and is on the low er half of the ran ge observed to- 
wards Galactic Hn regions dHonv et all 1200 ll) . The 12.7/11.2 
band ratio, on the contrary, is more similar to the ratios ob- 
served towards Galactic intermediate mass star-forming regions 
and Planetary Nebulae than towards Galactic Hn regions. This 
suggests the presence of more compact, smooth-edged PAHs to- 
wards N160 A2. However, a much more detailed study of the CH 
out-of -plane bending modes in Magellanic Cloud H n regions is 
warranted. 

The extended emission TIMMI2 slits si and s2 registered faint 
and extended [S iv] emission. Extended emission from other 
atomic lines was not detected with TIMMI2. Figure|6]shows the 
variation of the [S iv] line peak flux across these slits from south 
to north. This emission practically extends across the 20" of sep- 
aration between the YSO and Al, and is most likely associated 
with the central star cluster located between Al and A2. This 
cluster extends over an area ~ 9" x 8" and its brightest compo- 
nent, if it were a main sequence star, could be, as mentioned in 
Sect.[321 an 06.5 V. 

A more complete map of this extended [S iv] line emission 
is shown in Fig. Q This map has been created using the mul- 
tiple Spitzer-IRS SL observations that are available. It is clear 
that the [S iv] emission is stronger at the positions of the Al and 
A2 components, while it is practically null at the position of the 
YSO. It is also evident that there exists an extended component 
towards the west and south-west. 
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Table 4. Line fluxes in units of 10 -20 W cm" 2 . 



Line/ratio 


A (pm) 


N88A 


N160A1 


N160A2 


N160A-IR 


Instrument 


Aperture (") 


Notes 


[Am] 


7.0 


< 5 


0.7 ±0.1 
6±2 


5.2 ±0.5 
< 6 


1.2 ±0.4 

< 6 


Spitzer-IRS (SL) 
ISO-SWS 


a x a (c) 
14x20 


1, 2 


[ Ar in] 


9.0 


< 10 
10 + 5 


< 17 
8.1 ±0.5 
70 ± 10 


< 22 
32 ±5 
100 ± 15 


< 16 
4.1 ±0.2 
64 ± 10 


TIMMI2 

Spitzer-IRS (SL) 
ISO-SWS 


1.2x/3 <a) 
a x a <c) 
14x20 


1, 2 


[Siv] 


10.5 


34 + 2 

66 ±21 
79 ± 12 


57 ± 3 
56 ±4 

290 ± 45 


46 ± 4 
111 ±3 
122 ±61 
205 ± 30 


< 12 
14.8 ± 0.4 

170 ± 25 


1 1MM12 

Spitzer-IRS (SL) 
Spitzer-IRS (SH) 
ISO-SWS 


1.2 X /3 <a) 

a x a (c) 
4.7 x 11.3 
14x20 


1, 2 


[Nen] 


12.8 


< 16 

3.3 ± 0.5 

< 5 


< 16 
7.0 ±0.8 

112 ±35 


28 ±7 
45 ±3 
59 ± 10 
200 ± 40 


< 36 
11.4 ±0.8 

245 ± 45 


TIMMI2 

Spitzer-IRS (SL) 
Spitzer-IRS (SH) 
ISO-SWS 


1.2x/3 <a) 
a x a (c) 
4.7 x 11.3 
14x27 


1, 2 


[Nem] 


15.5 


71 ±3 
90 ±21 


160 ±70 
440 ± 70 


210 ± 37 (b) 
370 ± 65 


96 ±8 
575 ± 100 


Spitzer-IRS (SH) 
Spitzer-IRS (LL) 
ISO-SWS 


4.7 x 11.3 
10.2 x 10.2 
14x27 


1, 2 


[Sm] 


18.7 


22 ±1 
< 11 


110 ±20 
320 ± 50 


154 ± 22 
270 ± 50 


68 ±5 
215 ±40 


Spitzer-IRS (SH) 
Spitzer-IRS (LL) 
ISO-SWS 


4.7 x 11.3 
10.2 x 10.2 
14x27 


1, 2 


[Sm] 


33.5 


12 ±3 
< 25 


98 ± 14 
640 ± 320 


154 ± 22 
600 ± 160 


545 ± 140 


Spitzer-IRS (LL) 
Spitzer-IRS (LH) 
ISO-SWS 


10.2 x 10.2 
11.1 x 22.3 
20x33 


1,2 


[Nem] 


36.0 


8±3 
< 21 


70 ±35 


37 ± 12 
73 ±20 


72 ±20 


Spitzer-IRS (LH) 
ISO-SWS 


11.1 X22.3 
20x33 


1, 2 


[Nem]/[Nen] 


15.5/12.8 


21 ±3 
> 18 


4± 1 


4± 1 
1.8 ±0.5 


2.3 ± 0.6 


Spitzer-IRS (SH) 
ISO-SWS 


4.7 x 11.3 
14x27 




[Srv]/[Nen] 


10.5/12.8 


>2 

20 ±7 
> 16 


> 4 
8± 1 

2.6 ± 0.9 


1.6 ±0.4 

2± 1 
1.0 ±0.2 


1.3 ±0.1 
0.7 ± 0.2 


TIMMI2 

Spitzer-IRS (SL) 
Spitzer-IRS (SH) 
ISO-SWS 







(a) /3(N88A)=4'.'5, y3(N160 Al)=4'/7, j8(N160 A2)=8'.'0 and /?(N160 A-IR)=3'.'5. (b) This line flux might be uncertain, see Sect. H/2] 
(c) q(N160 Al)= 3 r7, a(N 160 A2)=7'/4 and a(N160 A-IR)=5'.'5. Notes-(l) The AOT02 ISO-SWS line fluxes of N160 Al have been published 
by I Vermeii e t al. (2002a). For N88 A and N160 A2, we use instead the line fluxes we have directly measured from the APT 01 ISO-SWS spe ctra 
centered on these sources. (2) The full ISO-SWS spectrum centered at the position of N160 A-IR has been published by Henning et al. ( 1993). 

The YSO N160A-IR The spectrum of N160A-IR, shown in 
the lower two panels of Fig. [5] is dominated by a deep silicate 
absorption, indicative of the high extinction (around 70 mag, 
see Sect. I4.31 > existent towards this object. The TIMMI2 spec- 
trum lacks emission from fine-structure lines. Once again, the 
overall agreement between the TIMMI2 and Spitzer observa- 
tions is striking, in particular in the case of the continuum emis- 
sion/absorption. On the other hand, it is clear that the dominant 
contribution to the ISO-SWS aperture is the ionized gas sur- 
rounding the YSO. Even for the smaller Spitzer-IRS SL aper- 
ture (5'.'5x5'.'5), a contribution of the surrounding ionized gas 
is present in the form of fine-structure lines. In the case of the 
TIMMI2 observations, one might expect slit losses of about 
50%. However, the TI MMI2 spectrum agree s remarkably well 
with the photometry of Epc htein et al.l (Il984l) obtained through 
a 10" aperture (see Fig. [8}. 

In addition to the silicate absorption bands at 9.7 and 18 pm, 
the Spitzer spectrum shows PAH emission at 6.2, 7.7 and 
11.2 pm, and even marginally at 8.6 pm (see panel d of Fig. [5]). 



Table 5. PAH fluxes in units of 10~ 20 WcnT 2 . 



Source 


6.0+6.2 


7.7 


8.6 


11.0+11.2 


12.7 


Instrument 


N88A 








0+15 (a) 




IRS-SH 


N160A1 


11 


22 


1.4 


12 


y 


IRS-SL 


N160A2 


77 


165 


18 


45 




IRS-SL 










5+39 


14 


IRS-SH 




y 


693 


y 


60+233 


y 


ISO-SWS 



y means that the band is present but it is very weak. <a) No 11.0 pm band 
present. Notes-Integration limits are [7.07,8.2] and [8.2, 8.83] for the 
7.7 and 8.6 pm PAH bands, respectively. The integration limits for the 
11.0+1 1.2 //m bands are [10.85, 11.6]. 



4.2. Physical conditions and abundance calculations 

The electron density can be determined from the ratio of two 
atomic fine-structure lines of the same ionic species when they 
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Table 6. Elemental abundances. 



Arcsecond 



Fig. 6. Spatial variation of the [S iv] 10.5 jim line peak across slit 
#1 (si) and slit #2 (s2) of TIMMI2. North is to the right. Three 
distinct zones (1,2 and 3) are indicated (see also the left panel 
of Fig. |2]i. 



YSO 



[SIV] 10.5yur 



10~ 6 W/m 2 /sr 

Fig. 7. [Siv] 10.5 fim line flux map of N160A. Its patchy ap- 
pearance is the result of the limited spatial coverage. The dashed 
lines are a few examples of the Spitzer-IRS SL apertures. The 
brightest stars of the central c luster (stars #24, #35 and #58 of 
iHevdari-Malaveri et all [2002b are also indicated [see electronic 
edition of the Journal for a color version of this figure] . 





N88A 


N160A1 


N160 A2 


Ne/H 


1.7 ±0.1 


7 + 3 


4.9 ± 0.8 


(xicr 5 ) 


0.94 ± 0.04 <a) 


6.7 ± 0.9 (a) 


5.73 ± 0.02 (a 




2.1 < c > 


6.33 (b) 


4.48 (b) 




2.2 (d) 






S/H 


1.2 + 0.2 


7 ± 1 


3.6 ±0.7 


(xlO" 6 ) 


1.2 ±0.3 (a) 


7.8 ± 0.5 (a) 


6.5 ±0.4 (a) 




3.2 (c) 


7±2 (b) 


5 ± 1 (b) 




2.2 (d) 






Ar/H 




5.2 ± 0.3 


8± 1 


(xl(T 7 ) 




14.8 ±0.4 (a) 


15.1 ±0.2 (a) 



References: (a) infrared : Vermei & van derHulsd d2002h: (b) optical: 
IHevdari-Malaveri & Testori jl986h : (c) optical: lKurt et alj dl999h : (d) op- 
tical: lTestoretal.U2003h . 



are e mitted from l evels with nearly the same excitation energy 
(e.g. iRubin et"ail 119941) . This ratio is sensitive to gas densi- 
ties approximately in between the critical densities of each line. 
High density gas can be probed using line ratios such as [S in] 
33.5/18.7 and [Nem] 36.0/15.5 /urn. On the other hand, ionic 
abundances can be determined from the measured strengths of 
the fi ne-structure and Hi recombination lines (e.g. lRubin et all 
[19881) . 

We will now determine the physical conditions of the ion- 
ized gas (in terms of electron density and elemental abundances) 
in N88 A, N160A1 and N160A2. The final abundances of Ne, 
S and At are listed in Table [6j where they are compared with 
previous optical and infrared determinations. 

4.2.1. N88A 

We use directly the line fluxes measured by Spitzer-IRS. De- 
reddening is not an issue because of the low extinction at these 
wavelengths. Note that the extinction in e.g. the near-infrared 
is only 2% of the attenuation suffered in the optical (which is 
around 1 .5) when ap plying the SMC extinction law obtained by 
iGordon et al.1 d2003l) . 

Assumin g an electron temperature of 14 000 K for 
N88A (e.g [Kurt et all [1991 IVermeii & van der Hulstt 12001 
iTestor et"all 1200 3). the observed [S in] 33.5/18.7 line ratio gives 
an electron density of 2275(+390, - 310) cm" 3 . This density 
is com parable to the one obtained by Vermeii & van der Hulsll 
d2002l) using ISO-SWS (2200 cm" 3 ), and to the rms electron 
density (3500 cm" 3 ) der ived from radio continuum observations 
dlndeb etouw et al. 2004 measured a flux density o f 91 mJy at 6 
cm) using Eq. A. 2. 5 of Panagia & Wal mslevl d 1 9781) and a source 
diameter of 3" (see Sect. 12. 11 1. The observed [Nem] 36.0/15.5 
ratio, on the other hand, indicates a density < 10 4 cm" 3 . 

In the calculation of ionic abundances, we use the H/3 
line flux (2.6 x 10" 11 ergs" 1 cm" 2 ) predicted directly from ra- 
dio continuum observations. I t is estimated using Eq. 10 of 
iMartfn-Hernandez et al.l d2005l) . When t his H/3 line flux is com- 
ared with the line flux measured by Hevdari-Mal averi et al.l 
1999), we get a visible extinction Ay ~ 1.8 which is in 
agreement with the values obtained by previous authors (see 
Sect. I3.lt . We also assume an electron density of 3500cm" 3 . 
This density is lower than the critical densities of the lines in- 
volved in the calculations and hence, is not a decisive factor. We 
obtain Ne + /H + = 1.5(0.2) x 10" 6 , Ne ++ /H + = 1.5(0.1) x 10" 5 , 
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S ++ /H + = 7.2(0.3) x 1(T 7 and S +3 /H + = 5(2) x 10" 7 , where the 
uncertainties are given within brackets. 

The total Ne abundance can be determined by adding the 
contributions of Ne + and Ne ++ . The final Ne/H we estimate 
agrees well with the optical me asurements. It is, however , highe r 
than the abundance obtained by Vermeii & van derHulstl (f2002). 
but these authors only give the contribution of Ne ++ since they 
could not determine the ionic abundance of Ne + . 

In the case of S, adding the contributions of S ++ and S 3+ 
only accounts for ~50% of the expected sulphur abundance, i.e. 
about 0.2 times the s olar S abundance (wh e re we use the solar 
S abundance given bv lAsplund et aUl2005l) . lTestor et ail d2003l) 
estimate a S ++ /H + abundance only 30% higher than the value we 
calculate and an insignificantly low S + /H + . It might then be in- 
ferred that 50% of the S abundance must be in the form of S +4 , 
which needs energies around 47 eV in order to be produced. In 
fact, the optical estimates, which use ionization correction fac- 
tors to account for unobserved ions, give higher S abundances. 



4.2.2. N160A1 

We use directly the line fluxes measured by Spitzer-IRS (SL 
and LL modes). For an e le ctron t emperature of 10000 K 
dHevdari-Malaveri & Testorl 119861: IVermeii & van der Hulstl, 
2002), the observed [Sm] 33.5/18.7 line ratio gives an elec- 
tron density of 780(+410, -250) cm' 3 , simi lar to the density 
estimated by IVermeii & van der Hulstl d2002l) using ISO-SWS 
(around 450 cm -3 ), but somewhat lower than the rms electron 
density of 2550 crrT 3 . This rms electron density is derived from 
radio continuum observations assuming a source diameter of 3" 1 
(see Sect. 12. Ill and usin g the 6 cm flux density measured by 
llndebetouw et alj d2004l) . 

In the calculation of ionic abundances, we use the H/3 
line flux predicted directly from radio continuum observations, 
and assume an electron density of 2550 cm -3 -this density is 
again lower than the critical densities of the lines involved in 
the calculations. This H/3 line flux (1.7 x 10" 11 ergs" 1 cm -2 ) 
is similar to the de- r eddene d H/3 line flux estimated by 
iHevdari-Malaveri et all d2002l) . For Ne we obtain Ne + /H + = 
6.1(0.7) x 10" 6 and Ne ++ /H + = 6(3) x 10" 5 . For S, S ++ /H + = 
6(1) x 10" 6 and S +3 /H + = 7.8(0.6) x 10" 7 . And finally, for Ar, 
Ar + /H + = 3.5(0.5) x 10" 8 and Ar ++ /H + = 4.9(0.3) x 10" 7 . 

The abundances of Ne, S and Ar are given in Table [6] The 
abundances of Ne and S agree well with the optical and infrared 
estimates. The Ar abundance we obt ain, however, is about 3 
times lower than the one estimated by IVermeii & van der Hulst 
(2002) using ISO-SWS. ISO presented a clear advantage over 
other instruments when estimating elemental abundances: fine- 
structure lines and strong Hi emission lines (e.g. Bra and Br/3 
at 4.05 and 2.62 yum, respectively) were observed using roughly 
the same aperture. Since in the estimate of ionic abundances we 
are assuming that the volume occupied by both ionic gases (i.e. 
H + and X + ', where X is a given element) is the same, aperture 
mismatches clearly have an effect on the results. This effect is 
more pronounced in the case of Ar + , probably because this ion 
has a lower ionization potential (27.6 eV) than others such as 
S ++ (34.8 eV) and Ne + (41.0 eV). The volume occupied by Ar + 
is thus largely extended, and we are likely underestimating the 
[Arn] 7.0/im line flux. 



4.2.3. N160A2 

As it has been done previously, we use directly the line fluxes 
measured by Spitzer-IRS. In particular, we use those measured 
in the SH mode, which cover a wider wavelength range and, 
when there is an overlap, agree very well with the line fluxes 
measured in th e SL mode. Assuming an electron temperature of 
10 000 K (e.g. IHevdari-Malaveri & Testorl IT986I) . the observed 
[Sm] 33.5/18.7 line ratio gives a much lower electron density 
(< 200 cm 4 ) than the rms electron density derived from the 
6 cm continuum observations, i .e. 1850 cm" 3 , est i mated by us- 
ing the 130 mJy measured by Indeb etouw et al.l d2004l) and a 
source diameter of 5'.'3 (see Sect. 12. U . We believe this is be- 
cause of the different apertures through which both [S in] lines 
are observed (see Table 0); we might perfectly be overestimat- 
ing the [S m] 33.5 fim line flux associated with A2. The observed 
[Nem] 36.0/15.5 indicates a density < 10 4 cm" 3 . 

In the calculation of ionic abundances, we use the H/3 line 
flux (4.6 x 10" 1 1 erg s" 1 cm" 2 ) predicted directly from radio con- 
tinuum ob servations. When compared w ith the H/3 line flux mea- 
sured by Hevd ari-Malaveri et al.l d2002l) . we get a visible extinc- 
tion around 2.5, slightly higher than the mean estimation given 
by these authors (~ 1.5, see Sect. 13. 2\ . We also assume an elec- 
tron density of 1850 cm" 3 . This density is again lower than the 
critical densities of the lines involved in the calculations. We ob- 
tain Ne + /H + = 1.9(0.3) x 10" 5 , Ne ++ /H + = 3.0(0.5) x 10" 5 , 
S ++ /H + = 3.0(0.4) x 10" 6 and S +3 /H + = 6(3) x 10" 7 . We 
can as well estimate the Ar + and Ar ++ ionic abundances us- 
ing the [Am] and [Arm] line fluxes measured by Spitzer-IRS 
in the SL mode. We obtain Ar + /H + = 1.0(0.1) x 10" 7 and 
Ar ++ /H + = 7(2) x 10" 7 . 

The abundances of Ne, S and Ar are given in Table|6] The Ne 
abundance agrees well with the optical and infrared estimates. In 
the case of S, adding the contributions of S ++ and S 3+ only ac- 
counts for ~50% of the exp ected sulphur abundance, similarly 
to what we find for N88 A. IHevdari-Malaveri & Testorl ([l986) 
estimate a S ++ /H + abundance only 15% higher than the value 
we calculate, and an extremely low S + /H + abundance (about 
3 x 10" 7 ). On the other hand, for this region the S +3 /H + is suf- 
ficiently low that one does not expect a significant contribution 
in the form of S +4 . Adding the contribution of S + , the total S/H 
abundance we obta in is ~ 4 x 10~ 6 , abo ut 0.3 times the sulphur 
solar abundance (Asplund et al., 2005), instead of the canoni- 
cal 0.5 value. This abun dance is in between the lower and up- 
per limits determined by IHevdari- Malaveri & Testorl dl986l) . It 
is, however, slightly lower tha n the S abundance obtained by 
IVermeii & van der Hulstl d2002l) . 

A larger discrepancy is found in the case of Ar. The Ar abun- 
dance we obtain is about 2 times lower than the one estimated by 
IVermeii & van der Hulstl d2002l) using ISO-SWS. As mentioned 
previously for N160A1 in Sect. 14.2.21 aperture mismatches are 
probably the cause of this difference. 

4.3. The spectral energy distribution of the massive YSO 
N160A-IR 

Table [7] compiles all the available photometry for N160A-IR, 
where no color corrections are applied. The IRAC photometry 
has been obtained using an aperture of 772. For MIPS, we use 
a 16" aperture (note that the pixel scale is 4" pi xel" 1 ) centered 
on the position of the YSO. The photometry of Epc htein et al.l 
d 19841) was obtained using apertures of 373 for the L and M 
bands, and 10" for the Nl, N2 and 7V3 bands. Far-infrared KAO 
flux densities are integrated over 1' and hence are taken as up- 
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Fig. 8. Spectral energy distribution of N 1 60 A-IR between 1 /vm and 1 . 3 m m. The Spitzer-IRS SL and LL (2 order) s pectra are plot- 
ted in orange. The leg end of the data is as f ollows: filled circles are data bv lHenning et alj 0998), open diamonds by Epc htein et al.l 
(1984) and crosses bv lWerner et"al] dl978h : filled squares correspond to the Spitzer-IRAC photometry, the filled star to the Spitzer- 
MIPS photometry, and filled triangles are the IRAS fluxes. The best fit obtained with the dusty code is plotted by a dashed line [see 
electronic edition of the Journal for a color version of this figure]. 



Table 7. Photometry of N 160 A-IR. 



A 


AA 


Flux density 


Band 


Reference 


Qua) 




(mJy) 






1.25 


0.30 


<0.08 


./ 


1 


1.65 


0.30 


<0.24 


H 


1 


2.15 


0.32 


1.1+0.2 


K' 


1 


3.56 


0.75 


50±2 


IRAC-3.6/mi 


2 


3.76 


0.70 


80±9 


L 


3 


4.52 


1.01 


190±3 


IRAC-4.5 urn 


2 


4.69 


0.68 


267±25 


M 


3 


5.73 


1.42 


547±6 


IRAC-5.8yum 


2 


7.91 


2.93 


852±22 


IRAC-8.0yum 


2 


8.38 


0.96 


845±272 


Nl 


3 


9.69 


1.50 


437±241 


N2 


3 


11.88 


6.53 


<123xl0 2 


IRAS-12//m 


4 


12.85 


1.38 


3027±558 


N3 


3 


18.06 


4.12 


6081 ±1680 


Q 


3 


23.80 


9.12 


<113xl0 3 


IRAS-25 fim 


4 


50 


20 


<230xl0 3 


KAO 


5 


60.85 


25.89 


<623xl0 3 


IRAS-60//m 


4 


71 


19 


22900 + 100 


MIPS-70/im 


2 


100 


50 


<300xl0 3 


KAO 


5 


102.44 


39.55 


<791xl0 3 


IRAS-60/im 


4 


1.3d3 




199+14 


SEST 


1 



Refere nces: (D lHenning et al.l dl998h : (2) this work; (3) Epchtei n et al 
' 19841) ; (4) IRAS Point Source Catalogue v.2.1; (5) IWerner etal 



1978). 



per limits. We also take IRAS values as upper limits. The beam 
width of SEST at 1.3 mm is 22". The photometry is plotted, to- 
gether with the Spitzer-IRS SL and LL (2 nd order) spectra, in 
Fig. [8] The LL spectrum (2 order), which is observed in a 10" 
aperture, was scaled down by a factor of 1 .6 to the level of the 
SL spectrum, observed in a smaller aperture of 5"5. The 5-20 pm 
photometry agrees well with the observed spectra. 

The SED of N 160 A-IR has been reproduced wi th the dust 
radiative transfer model dusty (llvezic et al.1 119991) as an at- 
tempt at deriving approximate quantitative information about 
this YSO. For dusty we just have to specify the normalized spec- 
trum of the central star, the dust composition mix and its radial 
distribution. The code computes the dust temperature radial dis- 
tribution and the emerging radiation field. The dust composi- 
tion can be chosen from a variety of grain types, but we decide d 
to use the standard composition (Weingart ner & DraineL l200lb . 
i.e. silicate and gra phite from iDraine & Lea (Il984l) . and amor- 
phous carbon from Hanner (1988). The relative proportion of 
each component is a free parameter. A more detailed spectrum 
than the one we have is necessary to investigate the presence 
of other types of grains. We can, however, assess the presence 
or absence of silicate, graphite and amorphous carbon from the 
SED plotted in Fig. [8] The depth of the absorption at 9.7 yum, for 
instance, is very sensitive to the relative proportion of silicate. In 
the absence of silicate, graphite will be responsible for most of 
the emission below 9 yum, wh ile amorphous carbon w ill mostly 
emit between 20 and 100 /mi dPlante & SauvageLl2002l) . 

We placed a blackbody with an effective temperature of 
50000 K in the center of a spherically symmetric dust cloud. 
The result is hardly dependent on the exact value of this tem- 
perature within a reasonable range. The geometry of the massive 
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YSO, on the other hand, is almost certainly not spherically sym- 
metric, but the present data are unable to constrain more detailed 
geo metries. We also assu me a standard MNR grain size distribu- 
tion dMathisetal .1119771) . 

The best model fit (dashed line in Fig. [8j uses a mixture with 
silicates (50%), graphite (30%) and amorphous carbon grains 
(20%), and a single power law of the form p(r) oc r _L0 for the ra- 
dial profile of the density in the dust cloud. The observations are 
reproduced with an Ay = 70 + 10 mag. The dust cloud extends 
from an inner radius of a few hundred (~ 200) AU, at which the 
dust temperature reaches 590 K, to an outer radius of about 1 pc 
(but note that the observable size of the source will depend on 
the selected wavelength). The observed SED is well reproduced, 
specially the silicate feature at 9.7 pm. This model gives an op- 
tical depth of the 9.7 pm silicate of 4.29, and a ratio of the fluxes 
at 9.7 and 18 /mi of 2.7 x 10 _1 . 

The absolute bolometric luminosity of the central source is 
the simplest parameter to derive from the model. We obtain a 
bolometric luminosity of 2 x 10 5 Lq- It i s derived by scaling the 
model SED to fit the 1.3 mm data point, and using a value for 
the distance to the LMC of 50 kpc (IStorm et alll2004 . This lu- 
minosity sugg ests a central s tar of spectra l type Q6 V an d abou t 
30 Mq (iMartins et all 120051) . The fit by iHenning et al] (l998) 
derived a higher bolometric luminosity (1.4 X 10 6 Lq), but their 
SED was scaled up to the ISO-SWS spectrum centered on the 
YSO, whose continuum flux level is about 10 times higher than 
that observed by TIMMI2 and Spitzer-IRS (see panel e of Fig . |5j ■ 



5. Conclusions 

We have presented MIR spectra of two Magellanic Cloud Hn 
regions: N88A in the SMC and N160A in the LMC. N88A 
is a compact (~ 3" in diameter) nebula, but N160A is com- 
prised of 3 bright components: two high-excitation blobs, Al 
and A2, and a young stellar object, N160A-IR. These spectra 
have been obtained with TIMMI2 on the ESO 3.6 m telescope 
and the Infrared Spectrograph (IRS) onboard the Spitzer Space 
Telescope. With these high-spatial resolution data, we have been 
able to separate for the first time, and from a spectroscopic point 
of view, the three components of N 160 A. 

The spectra show a wide variety in terms of continuum, lines 
and molecular bands strength. N88A, N160A1 and N160A2 
show a rising dust continuum, PAH bands and strong fine- 
structure lines. These spectra are typical of Hn regions. On the 
contrary, the YSO present in N160A shows a strong dust con- 
tinuum with deep silicate absorption. 

The influence of spatial resolution is investigated by means 
of comparing these TIMMI2 and Spitzer-IRS spectra with ISO- 
SWS observations. Observations of N88A indicate that the 
fluxes observed by ISO-SWS and Spitzer-IRS come, in spite of 
their large apertures, exclusively from the compact 3" source. 
We derived a [S m] electron density of 2275 crrT 3 , and elemen- 
tal abundances of ~ 1 .7 x 10~ 5 for Ne, and of ~ 1 .2 x 10~ 6 for S. 

Concerning N 1 60 A 1 and N 1 60 A2, aperture effects are read- 
ily observed. In this respect, it is worthwhile to note the large dif- 
ferences observed in, for example, the [Ne m]/[Ne n] line ratio, 
given that this ratio is generally used as a diagnostic for the effec- 
tive temperature of the central ionizing star(s). This ratio tends 
to decrease with aperture size, which reflects the complexity of 
the ionization structure of N160A. The variations seen in the 
[Nem]/[Ne n] line ratio, with the highest value found close to the 
HEBs, suggest that the two HEBs are ionized by very hot stars, 
while the extended emission is ionized by cooler stars probably 



located in the central cluster. This central cluster is situated be- 
tween Al and A2, and extends over several arcseconds. 

For Al, we derived a [S in] electron density of 780 cirT 3 and 
elemental abundances of ~ 7 x 10~ 5 for Ne, of ~ 7 x 10~ 6 for 
S and of ~ 5.2 x 10~ 7 for Ar. For A2, we derived elemental 
abundances of ~ 4.9 x 10~ 5 for Ne, of ~ 3.6 x 10~ 6 for S and 
of ~ 8.0 x 10~ 7 for Ar. In the case of A2, the derivation of the 
electron density using the [S in] lines suffers from aperture mis- 
matches. We also observed weak and extended [S iv] emission 
with TIMMI2, which is most likely related to the central stellar 
cluster located between Al and A2. The PAH emission bands in 
both Al and A2 show intensity ratios and band profiles typical 
of Magellanic Cloud H n regions. 

Finally, we modeled the spectral energy distribution of the 
YSO N160A-IR with the dust radiative transfer model dusty. 
The best model has a mixture of silicates (50%), graphite (30%) 
and amorphous carbon grains (20%) with a single power law de- 
pendence of the density with radius proportional to r _1 , where 
r is the radius. This corresponds to a visual magnitude of 
70+10 mag and a 9.7 pm optical silicate depth of 4.3. A bolo- 
metric luminosity of 2 x 10 5 Lq was derived, suggesting a central 
star of spectral type 06 V and about 30 Mq. Thi s luminosity is 
about 10 times lower than the previous estimate of lHenning et al.l 
(119981) . 
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